Two-level Hamiltonian of a superconducting quantum point contact 
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In a superconducting quantum point contact, dynamics of the superconducting phase is coupled 
to the transitions between the subgap states. We compute this coupling and derive the two-level 
Hamiltonian of the contact. 

One of the key features of superconducting quantum point contacts (SQPC) is the existence of subgap states 
(so-called Andreev states) whose energies depend on the phase difference across the contact : 
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j-jjrj. where A is the superconducting gap, t is the normal transparency of the contact. Each transversal mode propagating 
through the contact generates two such states (with opposite energies). Thus, at energy scales less than A, it is often 
convenient to describe the contact as a set of two-level systems. Further we assume for simplicity that we have only 
a single propagating mode (and therefore only two subgap levels). 

However, for describing dynamics of the contact at nonconstant a, the information on the energy spectrum (|l|) at 
| each value of a is not sufficient. Mathematically speaking, we need a connection on the bundle of Hilbert spaces over 
. the circle of possible values of a, and, more specifically, the projection of this connection onto the two-level subspace. 
0^ | Technically, it amounts to computing the "dynamic" matrix element 

g: /(°) = <qi£ii>, (2) 

where |0) and |1) are the two subgap states at a given value of a. This quantity defines coupling between the dynamics 
' of the superconducting phase and transitions in the two-dimensional subspace of subgap states. 

To illustrate this point, consider a single SQPC connected to a grain of finite capacity C. This system has been 
studied previously in the adiabatic approximation 13 and in the two-level approximation P,|l4|]. The two-level 



Hamiltonian for this system may be written as 

o . 

Q . where TV is the dimensionless potential of the grain and H (a) is a 2 x 2 matrix. At each a, the eigenvalues of H (a) 
must be given by ([!]). However, this does not fix the whole dependence on a. In our earlier work || we suggested 

■ ' Another candidate for Ho(a) might be 

Hi{*) = (V JLY (5) 
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Obviously, the latter choice of H (a) would lead to a physically different behaviour of the system (||), although H^(a) 
has the same eigenvalues as H (a). The choice (^) of H (a) appears natural from the point of view of perturbation 
theory in backscattering. Moreover, we claim that expression (|J) is exact for any value of r (within the model described 
below). We omitted the derivation of (|J) in our paper |)| and now fill this gap in the present note. 
We shall describe the one-channel SQPC by the one-dimensional Hamiltonian: 

H juii = Hsc + H scatt . (6) 
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where and \& are electron operators (L and R subscripts denote left- and right-movers, (3 =|, | is the spin index), 
A(ir) is the superconducting gap. We assume the following coordinate dependence of the gap: 

in other words, the absolute value of the gap A is constant across the contact, while the phase changes by a at x = 0. 

The scattering part of the Hamiltonian H scaU corresponds to elastic scattering at x = and is also quadratic in 
electron operators. We further disregard the nature of the scattering and describe it by means of a scattering matrix. 

The Hamiltonian is quadratic and may be diagonalized by operators linear in ^ and 'J. We shall further compute 
the operators 7J and 7j corresponding to the subgap states. These operators satisfy the Bogolyubov-de-Gennes 
equations (l^] 
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(Here we related and 7} using the spin-rotational invariance of the Hamiltonian) . 

Solving the equations (ft), we find that, away from x — 0, u^{x) and v fl (x) have the following form: 
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where 



If we define 



k = VA 2 - E 2 . 
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the Bogolyubov-de-Gennes equations take the form: 



e^u R + v R 


= 0, 
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= 0, 
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The scattering matrix at x = matches u^{x) and v^x) at x = ±0 
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The former equation describes scattering of electrons, the latter one — scattering of holes. The amplitudes a and 
b must satisfy the unitarity condition: \a\ 2 + \b\ 2 — 1 (\a\ 2 — t, \b\ 2 — r). The same scattering amplitudes in ( |T5| ) 
and ( |l6| ) only assume that the scattering is spin-independent. We also neglect the momentum dependence of the 
scattering amplitudes (the so-called "instant scattering" approximation). Finally, we remark that the phase of a has 
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the same meaning as the superconducting phase a. Therefore, without loss of generality, we may assume that a is 
real: a* = a = \/t. 

The condition that the homogeneous system of linear equations (p^)-(p^) has a solution, reduces to 



which immediately gives (0) for the energy E and 
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(We assumed, without loss of generality, that < a < ir and that E > 0). 

Solving the system (Hj)-dl9) enables us to compute the following commutator: 
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We shall use this commutator to compute the "dynamic" matrix element 

,dH 



as follows: 



r) 1 r) J-f I 

1(a) = <0| — 11) = =-{01-5-11) = ^ 

Oa hi — Eq oa 2E 



(19) 



(20) 



Since H is quadratic in fermionic operators, X(a) is just a number. A straightforward calculation (with normalized 
7-f) l\) S ives 



\X(a)\ = 



2E(a) 
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The phase of X(a) depends on the choice of phases of subgap state operators 7^ and 7J or, equivalently, on the 
relative phase of the two states |0) and |1). To fix this phase, we observe that our system is invariant under the 
combined time reversal and particle-hole symmetry. More specifically, this symmetry acts on operators as follows: 



L/3' 
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together with the complex conjugation of coefficients. The phase £ is adjusted depending on the phase of the backscat- 
tering amplitude b. (This is a modified version of the well-known symmetry of Bogolyubov-de-Gennes equations jlgl ). 
If we choose the Andreev states to be self-conjugate, then 



(0|A|0) = (l|ii) = o, /(a)=<0|^- 
oa oa oa 



1) is real. 
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(In analogy to the ordinary quantum mechanics with a real Hamiltonian: we may choose all eigenfunctions to be real, 
then the matrix elements of real operators will also be real.) 
With this choice of phases, from (pp|) ,(21), 
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The operator of the charge on the grain, projected onto the two-dimensional subbundle spanned by the states |0) 
and |1), takes in the basis {|0), |1)} the form 
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The Hamiltonian (||) in the basis {|0), |1)} takes the form 
H = 



E(a) 
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Now we want to perform the rotation to the "fixed" basis where the charge operator is simply 



Q = (27) 
oa 



This results in the Hamiltonian (||) with 



H (a) = U(a)( E( £ ) _E(a)) U ~ 1{a) > (28) 

where U(a) is a rotation: 

The angle of rotation (p(a) is given by 



H>{<*)= I (0\-^-\l)da= f da = I arctan (VF tan (30) 

Ca J 4 1-isin ( a /2) 2 V 2 



which upon substituting in (|28|)-(|29|) gives the result 

To summarize, we have replaced the multi-body superconducting system by the quantum-mechanical two-level 
Hamiltonian for the superconducting phase across the contact. This two- level approximation is appropriate whenever 
the system stays away from the upper continuum of excitations. The Hamiltonian (|J) loses its validity at points where 
the upper Andreev state touches the upper continuum (at a = 2nn). Once a particle reaches this point, it will pass 
to the vacant levels of the continuum instead of following the localized subgap levels fl6|| . This effect is important 
for many non-equilibrium problems, for example, those with constant voltage applied to the contact ra. For most 
equilibrium problems PJT^] and for some non-equilibrium setups 0JlJ] the Hamiltonian defined by Eqs (|3|), (||) may 
be used as the two-level approximation 
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[18] The Hamiltonian used in ref. jl4| is incorrect away from the transition region a ~ 7T. It is appropriate for the problem 
posed, but cannot be used for other problems involving transitions away from a w 7r. In such cases, the Hamiltonian 
should be used instead. 
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